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Abstract. A stochastic forest model of young and old age class trees is stud¬ 
ied. First, we prove existence, uniqueness and boundedness of global nonneg¬ 
ative solutions. Second, we investigate asymptotic behavior of solutions by 
giving a sufficient condition for sustainability of the forest. Under this condi¬ 
tion, we show existence of a Borel invariant measure. Third, we present several 
sufficient conditions for decline of the forest. Finally, we give some numerical 
examples. 


1. Introduction. In 1975, Antonovsky [1] introduced a mono-species forest model 
with two age classes of trees: 

i ^ = Pv - 7i'v)u - fu, 

\f^=fu-hv. 

Here, u and v denote the tree densities of young and old age classes, respectively. 
The parameters p, h and / are a reproduction rate, mortality of old trees, and aging 
rate of young trees, respectively; while 7 ( 1 ;) = a{v — 6 )^ -I- c is a mortality of young 
trees, which is allowed to depend on the old-tree density. In addition, o, b, c, p, / 
and h are assumed to be positive constants. 

It is not difficult to see that for any pair (uq,Vo) of nonnegative initial values Uq 
and To, the system (1) possesses a nonnegative and global solution. Furthermore, 
( 1 ) possesses nonnegative stationary solutions given by 
1.0 = ( 0 , 0 ) 

2. p- = {iih€[K,h*]) 
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3. P+ = {iihG{0,h*]) 

where /i* = , h* = . The stability and instability of these solutions 

depends strongly on the magnitude of the mortality h of old age class trees (see 
Table 1). 


h 

{0,K) 


(/i*, oo) 

o 

unstable 

stable 

glob, asymp. stable 

p+ 

stable 

stable 

- 

p- 

- 

unstable 

- 


Table 1. Stability and instability of stationary solutions of (1) 


On the basis of (1), Kuznetsov et al. [11] introduced a mathematical model of 
mono-species forest with two age classes which takes into account the seed produc¬ 
tion and dispersion. The third author studied that model with his colleagues (see, 
e.g., [3, 4, 5, 15] and [17, Chapter 11]). It is shown that h plays a crucial role in the 
asymptotic behavior of solutions. 

In the real world, the parameters in the model may be random variables due 
to unpredictability resulting from environmental, ecological and biological distur¬ 
bances. In principle, the deterministic forest model can not handle randomness. 
Investigating the role of fluctuation of parameters by using stochastic models should 
be an interesting problem in environmental and ecological sciences. 

As mentioned above, the asymptotic behavior of solutions to the deterministic 
forest model depends strongly on the magnitude of h. Therefore, in this paper 
we restrict ourselves to consider a stochastic forest model, where h is perturbed 
by (Gaussian) white noise. Since Gaussian white noises can be expressed as the 
generalized derivative of a Brownian motion, we make a substitution: 

h ^ h — adwt 

in (1), where {wt, t > 0} is a one-dimensional Brownian motion defined on a filtered 
complete probability space (fl, {J^t}t>o, P), and cr > 0 is the intensity of the white 
noise. Our stochastic forest model is then formulated by Ito stochastic differential 
equations in 

I dut = {pvt - [a{vt -hY +C + f]ut}dt, 

= ifut - hvt)dt + avtdwt. 

In this paper, we study the stochastic forest model (2). We prove existence of 
unique global solutions to (2) and then study their asymptotic behavior. On one 
hand, we present a sufficient condition for sustainability of the forest. Under this 
condition, we also prove existence of a non-trivial Borel invariant measure. On the 
other hand, we give several sufficient conditions for decline of the forest. The results 
are illustrated by a few numerical examples. 

To prove existence of non-trivial invariant measures to (2), a common method is 
to find four one-dimensional processes, namely and z;^, which satisfy two 

conditions: 

(i) Ut and Vt are bounded by these processes, i.e. u^{t) < ut < u^{t) and u^(t) < 
Vt < v'^(t) for 0 < t < oo 
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(ii) These four processes do not hit the boundaries in the sense that there exist 
e > 0 and M > 0 such that 

fe < u^{t) < v?{t) < M a.s.,0 < t < oo, 

|e < v^{t) < v‘^{t) < M a.s.,0 <t < oo 

However, this can not be done because 

liminfwt = 0 a.s. (see Theorem 2.3). 

>-oo 

To overcome this difficulty, we use the semigroup method presented in [6, 13]. 
First, we establish some estimates for the average of integrals of solutions (see 
(19) and (20)). Then, we construct a strongly continuous semigroup generated by 
solutions of (2). Using these estimates and a theorem in [13], we show that the 
semigroup enjoys an invariant measure. 

The organization of the present paper is as follows. Section 2 proves existence and 
boundedness of unique global nonnegative solutions to (2). Section 3 investigates 
sustainability of the forest and existence of a Borel invariant measure. To the 
contrary, Section 4 presents some sufficient conditions for decline of the forest. 
Finally, Section 5 gives some numerical examples. 


2. Global solutions. In this section, we prove existence of unique global nonneg¬ 
ative solutions to (2) and show boundedness of solutions. 

Put = {{u, v);u > 0,v > 0}. Denote by the closure of in K^. Put 


Mq = inf{M; pv — [7(11) + f]u < 0 for every v > 0}. 


Then, 


Mq = inf \u; , —r — u <0 for every u > o| 

L 7(n) -b / i 

pv 


= sup 


>0 a{v - 5)2 -b c-b / 

py^ab'^ + C + f 


( 3 ) 


y/a (y/ ab'^ + c + f - i/a5)2 -b c -b / 


For biological reasons, throughout this paper, initial values for (2) are taken from 

Let us first prove existence of unique global nonnegative solutions to (2). We use 
the following lemma. 


Lemma 2.1. Consider the one-dimensional stochastic differential equation: 


dxt = (02 — aiXt)dt -b axtdwt, 
Xt\t=o = xo > 0, 


( 4 ) 


where 01,02 and a are positive constants. Then, there exists a unique global solution 
of (4) such that 

(i) For 0 < 03 < 00, 


lim 

t—>-oo 


log(a3 -b Xt) 


= 0 


t 


a.s. 
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(ii) For 1 < 9 < oo and 0 < T < oo, there exists ag^T depending on 9 and T such 
that 

Ex( < ag^T, 0 < t <T. 

In addition, if 1 < 9 < 1 + then ag^x is independent of T 

(iii) limsup(_^ 3 o Xt = oo and liminft_>.oo a;* = 0 a.s. 

Since the proof of the lemma is quite easy, we omit it. 

Theorem 2.2. Let {uq, Vq) G Then, there exists a unique global solution (ut, Vt) 
of (2) such that {ut,Vt)\t=o = (uo,vo) and 

{ut,vt) G IR+ a.s., 0 < t < oo. 

In addition, if uq + vq > 0, then 

{ut,Vt) G a.s., 0 <t < oo. 

Proof. Since all the functions on the right-hand side of (2) are locally Lipschitz 
continuous, there is a unique local solution {ut,vt) defined on an interval [0,r), 
where r is a stopping time having the following property (see, e.g., [2, 7]). If 
P{t < oo} > 0, then r is an explosion time on {r < oo}, i.e. 

lim(|ut| -I- |z;t|) = oo a.s. on {t < oo}. 

t—*-T 

Therefore, it suffices to show that r = oo a.s. and that Ut > 0,Vt > 0 a.s. for 
0 < t < oo. To prove this, we use the method in [14, 16]. 

Consider the four cases of initial values. 

Case 1: uq = "yo = 0. This is a trivial case, since Ut = Vt = 0 a.s. for 0 < t < oo. 
Case 2; {uo,vo) G IR+. 

Let fco > 0 be a positive integer such that uq and vq lie in the interval [^,fco]- 
Denote 

ffk = [^,k] X [p/c], /c = 1,2,..., 

then Let us define a sequence of stopping times by 

Tk = ini {0 < t < T-{ut,Vt) i Hk} (5) 

with the convention inf 0 = oo. It is obvious that {t/c}^^^ is nondecreasing. Hence, 
there exists a limit Too of this sequence as k ^ oo: 

Too = lim Tj, < T a.s. 

fc—>-oo 

Let us prove that Too = oo a.s. Indeed, suppose the contrary, then there would 
exist T > 0 and 0 < e < 1 such that 

P{too <T}>e. (6) 

Consider a positive function V on which is defined by 

V{u, v) = u'^ + v^ — log u — log r;, u > 0, ?; > 0. 

The Ito formula gives 

dV{ut,Vt) = [LV]{ut,Vt)dt -I- {2av^ - a)dwt, 
where the infinitesimal operator L is given by 

[L-]iu,v) = + [pv - hiv) + f}u]^ + ifu-hv)^. 


( 7 ) 
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It is possibly seen that 

[LV]{u,v) =2{p + f)uv - 2[7 (t) + f]u^ + (cr^ - 2h)v'^ + 7(1;) 
pv fu ^ , (T^ 

u V 2 

In addition, there exist Mi > 0 (z = I, 2) such that 


[LV]{u,v) < MiV{u,v) + M 2 : (u, 


+ ■ 


Therefore, 

rtATk 


ntATk rt/\Tk rtATk 

/ dV{us,Vs) < / [MiV{us,Vs) + M 2 ]ds + / {2avl-a)dws. 

Jo Jo Jo 

Taking the expectations of the two sides of the latter inequality, we obtain that 


lEIJ(zZtATfc 1 ’^tATk ) 

rtATk 

<V{uo,vo) + M 2 {t^Tk) + M{E / V{us,Vs)ds 

Jo 

< [^(uo,'Co) + + Ml /" EI/('UsArfc,TsArJds, 0 < t < T. 

Jo 

The Gronwall inequality then provides that 

E^(wtArfc , WtAr J < [V {uq, Vq) + M2T]e^^^ 

< [V{uo, Vo) + M 2 T]e^^^, 0<t<T. ( 8 ) 

Hence, 

[V{uo,Vo) + M2T]e’^^'^ > EI/(zZTAr, , I^TAtJ 

^ E[l{T„„<T}^(MTArfc,'yTATj] 

= H'^{T^<T}ViUrk,Vrk)]- (9) 

On the other hand, (5) gives 

V{uTk,Vrk) > minj/c^ -logfc, -log 

= min - logfc, logfc + on {t^ < oo}. (10) 

Thanks to (6), (9) and (10), we observe that 

[H(uo,To) + M 2 T]e^^'^ > emin - log A:, log A: + 

Letting k —>■ oo, we arrive at a contradiction: 

oo > [H(uo,uo) +M2T]e^^^^ = oo. 

Thus, Too = T = 00 a.s. Furthermore, {ut,vt) S a.s. for 0 < t < oo. 

Case 3: uq > Q,vo = 0. 

Let A:o > 0 be a positive integer such that uq lies in [^, A:o]. Denote 

Hl = [l,k] 

and 

Tk = ini {O <t < T-Ut i Hi} , k = ko,ko + 1 ,... 
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Clearly, the sequence has a limit Too as /c —)■ oo: 

Too = lim Tfc < T a.s. 

k—yoo 

Let us first show that 


Vt > 0 a.s., 0 < t <Tk. (11) 

Indeed, due to the comparison theorem (see [8]), Vt > v^ a.s. for 0 < t < Tfc, where 
Vf is the solution of this equation: 

{ dVf = —hvfdt + av^dwt, 
u|t=o = Vo = 0. 

Obviously, Vj = 0. Thus, (11) follows. 

Let us now verify that 

there exists a > 0 such that EvtATk ^ ® for 0 < t < oo. In addition, for 
any T > 0, there exists /3 > 0 such that E7(v(ATfc) < ,5 for 0 < t < T. 

Indeed, by virtue of the first equation of (2), 

du r / N , .if Pv 1 

— = 7 z;) + / - U 

dt L7(v) + / J 

< biv) + f]{Mo - u), 0<t<T, 

where Mq is defined in (3). Solving this differential inequality, we obtain that 

Ut<Mo-{Mo-uo)e-db^i^')+n<i\ 0 < t < r. (13) 


( 12 ) 


Hence, 

Ut < maxjuo, Mq} a.s., 0 < t < r. (14) 

Using (14) and applying the comparison theorem for the second equation of (2), 
we observe that 

Vt < Vt a.s., 0 < t < T, (15) 

where Vt is the solution of the one-dimensional stochastic differential equation: 


dvt = [max{uo, Mq} — hvt]dt + avtdwt, 

Vt\t=0 = Vq. 


(16) 


Thanks to Lemma 2.1-(ii) and (15), it is easily seen that (12) holds true. 

Let us finally observe that 

T = oo and (ut,Vz)sK^ a.s., 0 < t < oo. 

In view of (II), it suffices to show that Too = oo. Indeed, suppose the contrary, 
then there would exist T > 0 and 0 < e < 1 such that 

P{too < T} > e. 

Consider a positive function H on (0, oo) defined by 

H{u) = V? — logit, u > 0. 

The Ito formula then gives 



dH{us) 



2ul'f{vb - 


PVs 

Us 


^0 


+ -f{vs)]ds 
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Taking the expectation of the two sides of this inequality and using (12) and (14), 
we observe that 

EH{usATk)< f {2puiax{uo,Mo}a + l3)ds 

Jo 

< {2pmax{uo, AIo}a + I3)T, 0<t<T. (17) 

By using (17) instead of (8), we repeat the same argument as in Case 2 to 
conclude that Too = oo a.s. 

Case 4: uq = 0,vo > 0. The proof for this case is similar to one for Case 3. 

By the above arguments, the proof of the theorem is complete. □ 

Let us now show boundedness for the density u of young age class trees and for 
moments of the density v of old age class trees. 

Theorem 2.3. Let {ut,vt) be a solution of (2) such that (ut,Vt)\t=o = {uq,vo) € 
Then, 

(i) supg<(.^oo^ where M* = max{uo,Mo} and Mq is defined in (3) 

(ii) limsup(_,.gc Ut < Mg a.s. 

(iii) For any 1 < 9 < 1 + ^, there exists ag > 0 such that 

limsupEv® < ag 

>-oo 

(iv) liminft^oo T't = 0 a.s. 

Proof. Clearly, (i) and (ii) follow from (13) and (14). Meanwhile, applying Lemma 
2.1 to the equation (16) and using the fact that Vt < Vt, (hi) and (iv) follow. □ 


3. Sustainability of forest. In this section, we present a sufficient condition for 
sustainability of the forest. Under this condition, we also show existence of a Borel 
invariant measure on for the system (2). 


3.1. Sustainability condition. Let us show that if the intensity of noise and 
mortality of old age class trees are small enough, then the forest is sustainable. 


Definition 3.1. The system (2) is said to be sustainable if for every initial value 
{uo,vo) G K+\{(0,0)}, the solution {ut,vt) satisfies 

limsupErij > 0 and limsupEiit > 0. 

t—^OO t—>oo 


Theorem 3.2. Assume that 


h < 


Pf 


ah"^ + c + f 


and < 2( 


PJ 


ah"^ + c + / 


- h) 


(18) 


and (uo,vo) G K^\{(0,0)}. Let (ut,vt) he the solution of (2) with (itt,Ut)|t=o = 
(uo,vo). Then, there exists e > 0 which is independent of{uo,vo) such that 

liminf- [ Vsds > e, liminf- [ v^ds > e a.s. (19) 


and 

..If* ..If* 

liminf - / Kusds > e, liminf - / Kvgds > e. (20) 

t-fOO t Jq t-fOO t Jq 

As a consequence, (2) is sustainable. 
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Proof. By virtue of (18), there exists a constant k > 0 such that 


< K < h, 


ab^+c+/ 

^2 ^ 2{p-Kh) 


Consider a function Q on R+\{(0,0)} defined by 

Q{u, v) = log(M + Kv). 

Theorem 2.2 and the ltd formula then provide that 


dQ{ut,Vt) = [LQ]{ut,Vt)dt 


anvt 


-dwt 


( 21 ) 


Ut + KVt 

where the operator L is defined in (7). After some simple calculations, we obtain 
that 


[LQ]{u,v) = 


{kJ — C — f)u + {p — Kh)v 


au{v — by 


U + KV 2{u + Kvy U + KV 


( 22 ) 


Thereby, by using the estimate (i) of Theorem 2.3, we observe that 
l0g(ut + KVt) 


= log{uo + KVo) - / 

Jo 

> log(uo + kvq) — aM* 

Jo 


* av^Us — 2abugVs 
Us + KVs 


ds - 


( k / - c - f - ab'^)us + {p- Kh)vs 


Us + KVs 
2 2 2 


-dws 


Us + KVs 


2{Us + KVs)^- 


g M* + KV, 


-ds 


Us + KVs 


-dws 


rt n 


+ 


/o 


(k/ - c - f - ab‘^)us + (|0 - Kh)vs 
Us + KVs 


2 2 2 


2{Us + KVs)"^- 

Let US show that there exists ei > 0 such that for all (u, v) S 


ds 


ds. 


( k / — c — f — ab^)u + {p — Kh)i 


2 2 2 


2 


(23) 


(24) 


U+KV 2{u + kv)‘^ 

Indeed, (24) is equivalent to that F{u,v) > 0 for all {u,v) G where 

F{u, v) =[2(k/ — c — f — ab^) — ei]u^ + 2 [k(k/ — c — f — a^) 

+ {p — kH) — Kei]vu + [2k{p — kK) — kJ'cj'^ — K'^eijv'^. 

Since it is easily seen that there exists a small ei > 0 such that the 

quadratic equation F{u, u) = 0 in the variable u has two non-positive solutions for 
every u > 0. Thus, F(n, u) > 0 for all {u,v) G 

Proof for (19). Due to (23), (24) and the fact that Vt < Vt, where Vt is the 
solution of (16), we have 


log(M* -I- KVt) aM* 


t 


M* 


-ds 


^ log(uo -I- kvq) £i 1 /■* crKV 
~ t 2 t Jq U + KV 


dws 


Put 


Nt = 


(TKVs 


Us + KVs 


-dWs. 


(25) 

(26) 
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Then, {-/Vt}o<t<oo is a real-valued continuous martingale vanishing at t = 0. Fur¬ 
thermore, {iVt}o<t<oo has a quadratic form given by 


<">• = i ^ " '■ 


/o (Us + KVs)^ 

The strong law of large numbers for martingale (see, e.g., [9, 12]) then gives 


lim —- = 0 
>-oo t 


a.s. 


(27) 


In the meantime, applying Lemma 2.1-(i) to the equation (16) and using Theorem 
2.3-(ii), we observe that 

< M* a.s., 0 < f < oo. 

Taking the limit as t —)■ oo of the two sides of (25), we hence obtain that 


lim inf - 
t—>-oo t 


*^,2 


aM*v 
M* 


ds > — 
“ 2 


a.s. 


Since 


< -jUs and 


< we conclude that 

j liminft^oo \ fg Vsds > a.s., 


lliminft^oo 


2aM* 

£l 

2a 


a.s. 


from which it follows (19). 

Proof for (20). Taking the expectation of the two sides of (25), we have 

-ds 


< lim inf 


t—yoo 


< lim inf 

t—yoo 


\og{M* + KVt)_ ^ aM* ^ 


M* 


t 

■ nEvt 


t Jg M* P KV 


aM* 1 


E,vds 


t K t 

here we used the estimate 

log(M* -b t) < M* X, 0 < X < oo. 

On account of Lemma 2.1, the solution v of (16) satisfies the estimate 

Ext < ai, 0 < t < oo, 

where ai is some positive constant. We thus have shown that 


lim inf - 

i—>-oo t 


Eucds > 


eiK 

2aM*' 


Meanwhile, taking the expectation of the two sides of the second equation of (2), 
it follows that 


lim inf — 

t^OO t 


Ertcds = lim inf 

t-^oo L t 


= lim inf — 
t—>oo t 


E^^t h f* ■ 
- 1 — / Evgds 


Kvsds > 


hnei 

2aM*' 


Therefore, (20) has been verified. As a consequence. 


limsupEut > 0 
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and 

limsupEiij > 0 . 

t—XX) 

This means that (2) is sustainable. We complete the proof. □ 

3.2. Existence of Borel invariant measure. Let us show existence of a Borel 
invariant measure of the Ito process (ut, Vt), which concentrates on some domain of 
under the assumptions in Theorem 3.2. 

Let P be the transition probability of {ut,vt)'. 

P{t,x,y,K) = P{('Ut,i;t) G K; {ut,Vt)\t=o = 

for 0 < t < oo, (x, y) G K+, and K G S(M^). It is well known that (see, e.g., [ 6 , 13]) 

(i) P{t,x,y,-) induces a strongly continuous semigroup {Pt}o<t<oo of operators 

on the space of bounded continuous functions: 

Ptfix,y)=[ f{i,r))P(t,x,y,d^dr)), / G 

(ii) P{t,x,y, •) induces a positive contraction [Pt] on the space of 

finite signed measures: 

[fiPt]{K) = [ P(t,x,y,K)y,{dxdy), y, G K G 

Jr2 

Definition 3.3. A Borel measure u on (i.e. a positive measure which is finite 
on any compact set of is said to be invariant with respect to {Pt}o<t<oo if for 
0 < t < oo and K G 

[uP*](iG) = u(iG). 

The following result is well known. 

Theorem 3.4 (Michael [13, Theorem 5.7]). Let X be a locally compact perfectly 
normal topological space. Let {Qt}o<t<oo be a strongly continuous semigroup on 
Cb(X) generated by a transition probability on {X, B{X)). If there exists a nonneg¬ 
ative function g in the space Co{X) of continuous functions with compact support 
such that 

Qtg(x)dt = oo, X G X, 

then there exists a Borel invariant measure for {Qt}o<t<oo- 
We are now ready to state our theorem. 

Theorem 3.5. Let (18) be satisfied. Then, {Pt}o<t<oo has a Borel invariant mea¬ 
sure whieh concentrates on some domain o/K^ n {{u,v);u < Mq}. 

Proof. To prove this theorem, we construct a function g G Co(M^) which satisfies 
the assumption in Theorem 3.4. 

On account of Theorem 3.2, we have 

If* e 

liminf- / lr„>Eii;fis > - > 0 a.s. 

t —>00 t Jq — 2 ^ 2 

Using Theorem 2.3-(iii) and the Holder inequality, for any 0 < 0 < there exists 
ne > 0 such that 

-<liminf- / E[l/„^> . lUsjds 

2 t—foo t Jq '.—21 
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< lim inf < 

t-loo t Jo 






<neliminf / P{iis > ^}ds\ 

t-s-oo Jo 2 y 


Thereby, there exists eo > 0 such that 


1 /■* e 

liminf- / Pjug > -}ds > Cq- 
t—VOO t Jq 2 


(28) 


On the other hand, by Theorem 2.3-(iii), there exists a > 0 such that Eut < a 
for 0 < t < oo. The Markov inequality then provides that 

2a 1 eo eo 

- « I - 


Hence, 


Put 


0<t<oo 


2 a 


inf P < Ut < — e ^ 1-• 


eo 


eo 


(29) 


Let us show that 


K = {(u,v);0 <u< M*, < V < —}. 

2 eo 


liminf - [ P{(ms,i;s) S K}ds > ^ > 0. 
t->-oo t Jo 2 


(30) 


Indeed, suppose the contrary, then there would exist an increasing sequence {t„} 
such that —>■ 00 as n —>■ 00 and 


1 


eo 


— ¥{{us,Vs) e K}ds < — 


n = 1,2,3,... 


Jo 


Theorem 2.3-(i) and (29) then give 


1 


I'n Jo 




ds 


JO 


t 


n Jo 


Vs < 


Vs < 




eo J t' 


n Jo 


e 2a 


ds 


Hds-^ 

eo J t 


P{(us,Us) e K}ds 


n Jo 


^0 ^ 0-1 
2 2 

Combining this and (28), we arrive at a contradiction: 

>1 — eo + eo = 1. 

Therefore, (30) holds true. 

Let us fix a nonnegative function g G (70(1]^+) such that 


ff(x,l/) = 


1, (x, y) € K, 

0, (ai,j/) G 
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where D K is some bounded open set of . In view of (30) , we have 

/ Psgix,y)ds= lg{^,r])P{s,x,y,d^dri)ds 
Jo Jo Jr^ 

> / / g{i,g)P{s,x,y,didy)ds 

Jo Jk 

t 


= / P{(lts, Us) G iir}c?S —>■ oo 

Jo 


as t —>■ oo. 


Thanks to Theorem 3.4, we conclude that there exists a Borel invariant measure v 
on for {Pt}o<t<oo such that u(iG) > 0. By Theorems 2.2-2.3, u concentrates on 


some domain of n {(m,u); u < Mq}. The proof is now complete. 


□ 


4. Decline of forest. In this section, we show decline of the forest when either the 
mortality h of old age class trees or the intensity a of noise is large. More precisely, 
if either 

Pf f{p + 2abM*) 


h > min 


or 


> 


c + / ’ ab"^ + c + f 
{p + c - h)"^ 


2 c 


then the forest falls into the decline. Here, M* is defined in Theorem 2.3-(i). 

Theorem 4.1. Let {ut,vt) be the solution of (2) with {ut,Vt)\t=o = (uo,vo) G 
K+\{(0,0)}. Assume that h > min | ^^a^'^+c+f ^ } • Then, as t ^ oo, Ut and vt 

converge to 0 in expectation, i.e. 

lim = lim Eut = 0. 

t—>-oo 

In particular, Ut and Vt converge to 0 in probability: 


lim F\ut > C} = lim Pjut > C| = 0, 

t—^OO t—^OO 


C > 0. 


Furthermore, 


lim P{(ut,ut) G H} = 0 

>-oo 


for any compact set A C 


^_l_. 


(31) 

(32) 

(33) 


Proof. Let us first prove that Ut and Vt converge to 0 in expectation. Consider the 
two cases of the mortality h. 

Case 1 : h > It follows from (2) that 


Eut < Mo + Jq [/oEms - (c + f)Eus]ds, 
Eut = Mo + /q [/E^s - /lEusJds. 


Since the functions wi and zu 2 defined by 

zui(X,Y) = pY-(c+f)X, W 2 {X,Y) = fX-hY 

are non-decreasing with respect to arguments Y and X, respectively, the comparison 
theorem applied to the latter system provides that 


Eut < Xt, 

Eut < yt, 


0 <t < oo, 
0 < t < oo, 


( 34 ) 
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where {xt,yt) is the positive solution to the linear system: 

^yt 

with {xt,yt)\t=o = (uo,vo). 

From this system and a fact that h > , we observe that 

=[p/-fe(c + /)]ai, < 0 , 

1 = [pf - h{c + f)]yt < 0. 

Hence, hxt + pyt and fxt + (c + f)yt are non-increasing as t increases. As a conse¬ 
quence, there exist two nonnegative constants /3i and (32 such that 

j limt^^(hxt + pyt) = Pi, 

\ limt_,.oo [fxt + {c + f)yt] = P 2 - 

It is then seen that 

f 1 ™ „ _ 0l{c+f)-P2P 

I Ilnlt-foo Xt — , 

1 1 ™ „ _ P'jh-Pif 

— ^c+f)-pf 

This means that h^c+f)-pf ') ^ stationary solution of (35). Substi¬ 

tuting this solution for {xt,yt) in (35), we obtain that 

f p{P 2 h - Pif) - (c -F f)[Pi{c + f)- P 2 P] = 0, 

\f[Pi{c + /) - P 2 P] - h{P 2 h - Pif)= 0. 

Solving this system of algebraic equations, we arrive at 

Pi = P 2 = 0 . 

Hence, 


lim Xt = lim yt = 0 . 

t—^oo t—^ca 

Combining (34) and (36), we conclude that 

lim Eu( = lim Eut = 0 . 

t—^OO t—¥00 

Case 2 : h > ■ From (2) and Theorem 2.3-(i), we have 

' /-i 

Eut <Uo+ [pEus -|- 2abE{usVs) — (ab^ + c + /)EMs](is 

<uo+ [{p + 2abM*)E,Vs — {ab"^ + c + f)^Us]ds, 

Jo 

Evt = uo -I- /q [fEus - hEvs]ds. 

Using the same argument as in Case 1, we conclude that 

lim Eut = lim Eiij = 0 . 

t—>-00 t—XX) 

Let us now verify (32) and (33). For any 0 < Ci < C 2 ,0 < di < ^ 2 , 
E{{ut,Vt) G [ci,C 2 ] X [di,d 2 ]} < ]P{wt > Cl}, 


(36) 


P{ut > Cij < —Eut, Pjut > Cl} < —Eut. 

Cl Cl 


and 
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This together with (31) derives (32) and (33). It completes the proof of the theorem. 

□ 

Under somewhat stronger assumptions than those of Theorem 4.1, we can show 
almost sure convergence of Ut and Vt to 0. Consider two functions Fi and F 2 defined 
by 


Fiix) =fx^ + 2f{a^ + h-c- f)x^ + [(c + / - hf 



- 2 p/ - 2 (c -k f)a'^]x^ + 2p{c + f - h)x + 

(37) 

and 

^ 2 ( 2 ;) = fx"^ - {c + f + h)x + p. 

(38) 

Assume that either 



inf Fi{x) < 0 , 

(39) 

or 

2p c + f , , , , , , 

^5 -- < —-— and there exists A such that 

-I- 2 a 7 

< A < Ui(A) > 0 and F 2 (A) < 0 

cr^ - 1 - 2 a 7 

(40) 


holds true. Then, the following theorem shows such convergence. 


Theorem 4.2. Let {ut,vt) be the solution of (2) with {ut,Vt)\t=o = (uo,vo) S 
K+\{(0, 0)}. Under (39) or (40), limt_,.oo Ut = limt_>oo "Ct = 0 a.s. 

Proof. We again use the function Q defined by Q(u, v) = log(u + kv) as in the proof 
of Theorem 3.2, where k is a positive constant that will be fixed below. 

Let us first show that under (39) or (40), there exists a small e > 0 such that 

[LQ]{u,v) ^ foi' 3'b {u,v) S 

where [LQ] is defined in (22). Indeed, it is easily seen that a sufficient condition for 
this (in fact, it is also a necessary condition) is that there exists e > 0 such that 

F{u, v) =[2(c + / — Kf) — e]u^ — 2[/«(k/ — c — f)+p — Kh + Ke]vu 

(41) 

+ + 2 k(k/i — p) — n^efv^ > 0 for all (u, v) G 

If (39) takes place, choose k such that 

c f 

0 < K < — j — and Fi{k) < 0. 

It is then easily seen that there exists a small e > 0 such that the quadratic equation 
F(u,v) = 0 in u has a non-positive discriminant for all u > 0. This implies (41). In 
the meantime, if (40) takes place, choose k = A in (40). Similarly, it is seen that 
there exists e > 0 such that the equation F{u, v) = 0 has non-positive two solutions 
for all u > 0. This also derives (41). 

Let us now verify that 

lim Q{ut,vt) = —00 a.s. 

t—¥oo 

It follows from (21) that 

111 1 (7 KiV 

-Q{ut,vt) =-Q{uo,vo) + - / [LQ]{us,Vs)ds + - / - —^dws 

t t t Jn t Jn Us + KVs 
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1 If 

< -Q{uo,vo) - e + T / 

I t Jo 

On account of (26) and (27), 


auvs 


0 “s 


Us + KVs 


-dw.s 


0 <t < oo. 


(42) 


lim - 

t—^OC) t 


(TKVs 


-dws = 0 


a.s. 


/o Us + KVs 

Hence, taking the limit as t —>■ oo of both the hand sides of (42), we observe that 


limsup yQ(ut,Wt) < 
t—¥00 i ^ 


a.s. 


This implies that limt_>oo Q{ut,vt) = —oo a.s. Thus, 


lim ut = lim uj = 0 a.s. 

t—^OO t—>-oo 

The proof is complete. □ 


Remark 1. It is possibly seen that if Ti(l) < 0, then (39) takes place. After some 
simple calculations on the inequality T’i(l) < 0, we arrive at this condition; 


^2 ^ jp + c-h)^ 

2 c 

According to Theorem 4.2, we therefore conclude that a noise with large intensity 
causes decline of the forest. 


5. Numerical examples. Let us exhibit some numerical examples for sustainabil¬ 
ity of the forest and possibility of decline. For the computations, we used a scheme 
of order 1.5 (see, e.g., [10]). 

5.1. Sustainability of forest. In the system (2), set a = 2, 6 = 1, c = 2.5,/ = 
4, /i = 1, p = 5, (T = 0.5, and take initial value (uq, vq) = (2,1). 

Figure 1 gives sample trajectories of u and v in the phase space and in time. 
Figure 2 plots points {ut,vt) of 10^ sample trajectories of (u,v) at time T = 
1000 . 

By computing 10^ sample trajectories of {u,v), Figure 3 shows a graph of the 
expectation of tree densities of young and old age classes. 

Figure 4 gives a sample trajectory of two processes I and J defined by 

I{t) = \ [ Usds, J{t) = ^ f Vsds. 

Jo ^ Jo 

Figure 5 demonstrates a trajectory of two probability functions R and S defined 
by 

f R{t) = P{(ut, vt) e A; {uo, Vo) = (2,1)}, 

[Sit) = F{iut,vt) e A; iuo,vo) = (3,4)}, 

along t G [50,100], where A = [0.5,30] x [0.1,20]. These functions are calculated 
on the basis of 2000 sample trajectories of {ut,vt) corresponding to each of the two 
initial values. 
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Figure 1. Sample trajectories of Ut and Vt of (2) with param¬ 
eters: a = 2,6 = l,c = 2.5,/ = 4 ,6 = l,p = 5, cr = 0.5 and 
initial value (uo,vo) = (2,1). The left figure illustrates a sample 
trajectory of (ut,Vt) in the phase space; the right figure illustrates 
sample trajectories of Ut and Vt along t e [0,100]. 



Figure 2. Distribution of (ut,vt) of (2) at t = 10^. The parame¬ 
ters and initial value are taken as in the legend of Fig. 1. 
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5.2. Decline of forest. First, set a = 3, 6 = 4, c = 5,/ = 6,/i = 2,/9 = 7, cr = 4 
and take (uq, vq) = (4,3). Figure 6 gives sample trajectories of u and v in the phase 
space and in time. 

Second, set a = 3,6 = 4, c = 5,/ = 6,h = 3.82, p = 7, a = 0.25 and take 
{uq, Vq) = (4, 3). By computing 5 x 10^ sample trajectories of (m, u). Figure 7 shows 
a graph of expectation of tree densities of young and old age classes. 



Figure 3. Graphs of Eu and Ev along t G [0, 20]. The parameters 
and initial value are taken as in the legend of Fig. 1. 



Figure 4. Sample trajectory of two processes I and J defined 
by /(t) = j J^Usds and J(t) = j f*Vsds along t G [0,100]. The 
parameters and initial value are taken as in the legend of Fig. 1. 





















Probability 
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Figure 5. Graph of probability functions R and S defined by 
R{t) = P{(nt,nt) G A]{uo,vo) = (2,1)} and S(t) = F{{ut,vt) G 
A-,{uo,vo) = (3,4)1 along t G [50,100], where A = [0.5,30] x 
[0.1,20] and the parameters of (2) are taken as in the legend of 
Fig. 1. These functions are calculated on the basis of 2000 sample 
trajectories of {ut,vt) corresponding to each initial value. 





Figure 6. Decline of forest under the effect of noise with large 
intensity a. Here, a = 3, 6 = 4, c = 5,/ = 6, h = 2,p = 7, ct = 4 and 
initial value (uq, vq) = (4, 3). The left figure is a sample trajectory 
of {ut, Vt) in the phase space; the right figure is a sample trajectory 
of u and v along t G [0,1]. 
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Figure 7. Decline of forest when the mortality h of old trees is 
large. Here, a = 3, b = 4, c = 5, f = 6, h = 3.82, p = 7,a = 0.25 
and initial value (uo,vo) = (4,3). The figure gives a graph of Eu 
and Eu along t G [0,10]. 
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